Wannier-Stark ladder in the linear absorption of a random system with scale-free 

disorder 
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We study numerically the linear optical response of a quasiparticle moving on a one-dimensional 
disordered lattice in the presence of a linear bias. The random site potential is assumed to be long- 
range-correlated with a power- law spectral density S(k) ~ l/k a , a > 0. This type of correlations 
results in a phase of extended states at the band center, provided a is larger than a critical value 
a c [F. A. B. F. de Moura and M. L. Lyra, Phys. Rev. Lett. 81, 3735 (1998)]. The width of the 
delocalized phase can be tested by applying an external electric field: Bloch-like oscillations of a 
quasiparticle wave packet are governed by the two mobility edges, playing now the role of band 
edges [F. Dommguez-Adame et al, Phys. Rev. Lett. 91, 197402 (2003)]. We demonstrate that the 
frequency-domain counterpart of these oscillations, the so-called Wannier-Stark ladder, also arises 
in this system. When the phase of extended states emerges in the system, this ladder turns out 
to be a comb of doublets, for some range of disorder strength and bias. Linear optical absorption 
provides a tool to detect this level structure. 

PACS numbers: 78.30.Ly; 71.30.+h; 36.20.Kd 



INTRODUCTION 



Since 1998^ there exists an increasing interest in 
studying the localization properties of the quasipar- 
ticle wave functions in one-dimensional (ID) disor- 
dered systems with a long-range-correlated site poten- 
tial landscapei 2 . 3 . 4 . 5 . 6 . 7 . 8 . 9 i 10 . n i 12 i 13 i 14 i 15 i 16 i 17 i 18 i 19 . 20 i 21 i 22 
Random sequences, having a power-law spectral den- 
sity S(k) ~ l/k a with a > 0, result in a phase of ex- 
tended states at the band center, provided a is larger 
than a critical value a c m^ This finding contradicts the 
widely admitted conclusion of the one-parameter scal- 
ing theory of localizations 3 , that all states of noninteract- 
ing quasiparticles in one and two dimension s with time 
reversal symmetry are localized (see Refs. I24l25l26l27l 
for an overview). As a matter of fact, a great deal of 
work is being devoted to put the correlation-induced low- 
dimensional localization-delocalization transition (LDT) 
on solid grounds. 

It is worthwhile to notice that long-range correlations 
with spectral density of the form l/k a are widely pre- 
sented in nature, both in vitro and in vivo (see, e.g., 
Refs. and I2H and references therein). This type of 
correlation gives rise to the " fractal geometry of nature" 
introduced by Mandelbrot. 30 Importantly, in 1992 it was 
conjectured that the long-range power-law correlations 
exist in nucleotide sequences in dnAj 3 ^ 3 ^ 3 ^ 3 ^ This con- 
jecture opened prospectives to quantify nature in vivo 
with critical exponents^ 5 - It was argued also that long- 
range correlations in DNA sequences can explain the 
long-distance charge transport in these systems (ASiiiiSfliSi 
The l/fc a -law has also its trace in energy level statis- 
tics^ All said above unambiguously testifies that study- 



ing the properties of disordered systems with long-range- 
correlated disorder is an attention grabbing task. 

It is known since seminal papers by Bloch- 7 and 
Zener— that an electron moving in a periodic potential 
and subjected to a uniform electric field is localized due 
to the Bragg reflection. It performs a periodic motion, 
known as Bloch oscillations ^SMS which is characterized 
by an angular frequency lob — eFd/Ti and a spatial ex- 
tension Lb = W/(eF), where — e is the electron charge, 
F is the applied electric field strength, d denotes the 
spatial period of the potential, and W stands for the 
band width. The Bloch oscillations were observed for 
the first time as oscillations of electronic wave-packets 
in semiconductor superlattices^i^iM 4 * 4 ^ 6 . (see for an 
overview Ref. 147]) , and later on as a periodic motion of 
ensembles of ultracold atoms^f^ 9 . and Bose-Einstein con- 
densates 50 in tilted optical lattices. The multiplicity of 
observable physical phenomena related to electron Bloch 
oscillations increases even more when a semiconductor 
superlattice is subjected to joint, perpendicular or tilted, 
electric and magnetic fieldsi 51 . 52 . 53 . 54 . 55 

Recently, it was demonstrated that ID disordered sys- 
tems with the l/k a spectral density support Bloch-like 
oscillations of quasiparticles i& It was also shown that 
these oscillations provide a tool to measure the energy 
width of the delocalized phase arising at a > a c . More 
specifically, the two mobility edges, which separate the 
phase of extended states from the two phases of local- 
ized ones, were found to play the role of effective band 
edges, i.e., it is the width A of the mobility band who 
determines the spatial extension L* B = A/ (eF) of Bloch 
oscillations. Therefore, this width can be measured in 
biased disordered lattices. 
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In this work we report further progress along this 
line by considering the frequency-domain counterpart of 
the Bloch oscillations, the so-called Wannier-Stark lad- 
der (WSL)£& (see Refs. H3 and |H for brief histori- 
cal surveys). The WSL is characterized by a series of 
equidistant quasistationary levels separated by an energy 
U — Tiljb — eFd. The progress in semiconductor growth 
techniques made it possible to firmly establish the exis- 
tence of ladder level structures in semiconductor super- 
latticea 59 i 60 i 61 i 62 as well as in <5-doped super lattices ■ 63 1 9 4 
There exists evidence that moderate uncorrelated disor- 
der does not destroy the required phase coherence to see 
the WSL in continuous (many band) models^ We focus 
on the one-band approximation and demonstrate that in 
our system the WSL also arises, in spite of the underlying 
randomness. Strong correlations in the disorder facilitate 
the observation of the WSL. At a > a c , when the phase 
of extended states emerges at the center of the band, the 
WSL may present a comb of doublets, reflecting the dou- 
blet energy structure of the unbiased systemi£& To work 
out the problem, we numerically calculate the absorption 
spectrum varying the correlation exponent a, the disor- 
der strength, and the magnitude U of the bias. 

The outline of the paper is as follows. In the next 
section, we present our model which is based on a tight- 
binding Hamiltonian of a quasiparticle, moving in a long- 
range-correlated potential landscape and subjected to a 
linear bias. In Sec. IIIII we recall the basic physics of the 
linear absorption spectrum in the absence of both bias 
and correlations in disorder. The central parts of the 
paper are Sec. IIVI and Sec. E] where the results of nu- 
merical simulations of the absorption spectrum profile in 
disorder-correlated systems are shown. We begin with 
a brief discussion of the absorption spectrum behavior 
upon increasing the disorder correlations in bias-free sys- 
tems rSeclIVjl. proceeding in Sec.EJonto the absorption 
in biased lattices. We discuss in detail its dependence on 
the driving parameters of the model (bias magnitude, dis- 
order strength and correlation exponent a) and provide 
an evidence of that the correlations in disorder facilitate 
the occurrence of the WSL. Finally, Sec. IVII concludes 
the paper. 



II. MODEL 

We consider a biased tight-binding model with diago- 
nal disorder on an otherwise regular ID open lattice of 
spacing unity and N sites (N is assumed to be even) . Wc 
assign to each lattice site two levels (ground and excited 
states) with transition energy £ n and consider optical 
transitions between them. The model Hamiltonian is 



JY 



w = E 



n=l 
N-l 



\n)(n\ 



£ (|n)(n + l| + |n + l)(n|) 



(1) 



Here, \n) denotes the state in which the n-th site is ex- 
cited, whereas all the other sites are in the ground state. 
The energy of the state \n), 8 n = £ + e n , is assumed to 
have a stochastic part e n generated according toi 
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e n = °C a ^-^ cos — 

k=l V 



(2) 



where C a — \/2( X^fcfi ^ 1 2 is the normalization 
constant and <pi, . . . , 4>m/2 are N/2 uncorrelated random 
phases uniformly distributed within the interval [0, 2ir]. 
The distribution has zero mean (e n ) = and standard 
deviation (e 2 ) 1 / 2 = cr, where (...) indicates averaging- 
over realizations of random phases <pk- The quantity a 
will be referred to as magnitude of disorder. The stochas- 
tic sequence is characterized by a correlation function 
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which is long-ranged, except for the particular value of 



the exponent a — 0, when (e n e m ) = a 2 Sn 



being 



the Kronecker S. The site energies in this case are uncor- 
related. Correlations, however, arise as soon as a =/= 0. 
Thus, for a < 1, they are power-law-like, i.e., the cor- 
relation function JSJ decays as \n — m| Q_1 . If a ^> 1, 
the term with k = 1 in the series © is dominant, and 
(e n e m ) = (l/2)a 2 C 2 cos[27r(n— m)/N], implying full cor- 
relation of the on-site energies. 

The term —U(n — N/2) in the Hamiltonian de- 
scribes the linear bias. We will not necessarily relate 
it to the presence of an external uniform electric field, 
as in the case of an electron moving in the conduction 
bandS The bias can also be an intrinsic property of the 
system as it takes place in dendritic species (see, e.g., 
Ref . l68l and references therein). In this case, the Hamil- 
tonian models a ID Frenkel exciton in a disordered 
lattice with energetic bias. Finally, the intersite transfer 
integrals in JQ) are restricted to nearest-neighbors, and 
it is set to —1 over the entire lattice. Also, we set £ = 
hereafter without loss of generality. 

As we already mentioned in the Introduction, in the 
absence of bias (U = 0), the above model supports a 
phase of extended states at the center of the band, pro- 
vided the correlation exponent a is larger than a crit- 
ical value a c . At a < a c all the states are localized, 
which implies that the model under consideration under- 
goes an LDT with respect to the correlation exponent a. 
In Ref. U where the above model of disorder was intro- 
duced, the disorder magnitude was set to a = 1, and the 
critical value a c was found to be q c = 2. It may seem 
that a c depends on a. However, it was demonstrated 
later on that a c = 2 is the universal critical value for the 
LDT to occur in this model, independently of 

Another peculiarity of this model, having a direct rela- 
tionship to the specific form of the random potential 
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as well as to its localization properties, is that the ab- 
sorption spectrum at a > a c = 2, i.e., in the presence 
of the phase of extended states, reveals a double-peaked 
structure (see Ref . lo6l and Sec. lIV|) . 



III. ABSORPTION SPECTRUM 

The quantity subjected to calculation throughout this 
paper will be the absorption spectrum defined as 



A(E) 




8{E-E y ) 



(4) 



where the E v and Vw are the eigenenergies and eigen- 
functions, respectively, obtained after diagonalization of 

the Hamiltonian l|TJl. The quantity F v = (j")^^ tp U nj 

is the dimensionless oscillator strength of the v-th state. 

In order to gain insight into the effects of long-range 
correlations in disorder and the bias on the absorption, 
we first recall the basic features of the absorption spec- 
trum in the absence of both (a — and U = 0). Un- 
correlated diagonal disorder results in localization of all 
the states^ and in the appearance of Lifshits tails in the 
density of states (DOS), outside the bare quasi-particle 
band which ranges from E = —2 to E = 2. Since we set 
negative sign of the hopping integral, the low-energy part 
of the DOS (with E around —2) is of importance for the 
linear optical absorption. The majority of states lying in 
this region are localized at different (weakly overlapped) 
segments. Some of them are bell-like, i.e., without nodes 
within localization region, while the other and higher 
(band) states resemble standing waves with nodesiS The 
bell-like states dominate the optical absorption because 
they accumulate oscillator strengths large compared to 
those of the other states. They result from localization 
on fluctuations of the site potential, which have a well- 
like envelop^ The typical size of a such potential wells, 
N* , determines the extension of the bell-like states, while 
the depth a/y/N* governs the width of the absorption 
spectrum (see below). Notice that the potential depth is 
y/W* times as small as the bare magnitude a. This effect 
is known as the exchange (or motional) narrowing: fluc- 
tuations of a stochastic site potential of alternating signs 
are averaged out by a quasiparticle rapidly moving (due 
to a large exchange interaction J 3> c) within a region 
of size _/y*, 69 i 71 i 72 i 73 i 74 

The absorption spectrum is peaked slightly below the 
bare band edge E = —2 and represents an inhomoge- 
neously broadened line with a Gaussian-like red and a 
Lorentzian-like blue tail (see, e.g., Ref. Fl|) . The local- 
ization size of the optically dominant states, N* , and the 
full width at half maximum (FWHM) of the absorption, 
a*, are estimated as— 



N* = 



3^ 
er 



2/3 




FIG. 1: Absorption spectra of an ensemble of unbiased chains 
(U = 0) with N — 500 sites calculated for two values of the 
correlation exponent a, shown in the plot. In both cases, the 
magnitude of disorder is a = 1. Each curve were obtained 
after averaging over 3 x 10 4 realizations of disorder. 
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To conclude this section, note that both short-range 
and power-law long-range correlations in disorder result 
in decreasing the localization length of the band edge 
state d 22 ! 74 and, subsequently, increasing the absorption 
line width i 66 ! 74 In contrast, this is not the case for cor- 
relations which are stronger than power-law-like. An ex- 
ample is a sequence © at a ^> 1 (see the next section). 



IV. UNBIASED SYSTEM 

We first discuss the absorption spectrum of the unbi- 
ased system (U — 0), aiming to separate the effects of 
bias from those related to the disordered nature of the 
model. Figure H represents the absorption spectra calcu- 
lated for two values of the correlation exponent, a = 1 
and a = 4. By convention, we will refer to these two 
cases as weakly and strongly long-range correlated disor- 
der, respectively. The results were obtained by numeri- 
cally diagonalizing the Hamiltonian for chains of size 
N = 500 with open boundary conditions. The disorder 
magnitude was set to a = 1. Averaging over 3 x 10 4 re- 
alizations of the disorder were applied in Eq. Q for each 
value of a. 

From Fig. ^ we observe that in the case of weak cor- 
relations in disorder (a = 1), the absorption spectrum 
consists of a single inhomogeneously broadened asym- 
metric line. Its red and blue tails can be fitted by a 
Gaussian and Loretzian, respectively. The spectrum is 
peaked slightly below the low-energy band edge E = — 2 
and has the FWHM a* w 2.2 which is approximately 
three times larger than that in the limit of uncorrelated 
disorder (a = 0), Eq. Ij5b|l . This tendency is in full accor- 
dance with that we mentioned in the preceding section. 

The absorption line shape, however, changes dramat- 
ically when disorder is strongly long-range correlated 
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(a = 4). The higher-energy peak in the absorption spec- 
trum starts to build up when the correlation exponent a 
exceeds the critical value a c = 2 or, in other words, when 
a phase of the extended states emerges in the center of 
the band. 66 Its appearance implies that states deep inside 
the band gain large oscillator strengths. This is in con- 
trast to the case of uncorrelated disorder, where higher 
states have a vanishingly small oscillator strengthiiiS 

A simple explanation of this anomaly is based on the 
fact that for sufficiently large values of a, the first term in 
the series © is dominant, while the others are consider- 
ably smaller. Consequently, the site potential for a given 
realization is cosine- like (harmonic with k = 1), per- 
turbed by a colored noise (harmonics with k > 2). Then, 
the whole lattice can be represented as a two weakly- 
coupled sublattices with different site energy (see Ref. |S£j 
for more details). Optical transitions to the band edge 
states of these two sublattices give rise to a double-peaked 
structure of the absorption spectrum. Remarkably, the 
higher-energy peak monitors the upper mobility edge of 
the delocalized phased 



V. BIASED SYSTEM 
A. Qualitative picture 

In disorder-free systems, switching the bias on re- 
sults in dynamical (Bloch) localization of all the states^ 
within the localization size Lb = 4/f/, where 4 is the 
band width in dimcnsionlcss units. The Bloch localiza- 
tion is accompanied by the subsequent reorganization 
of the energy spectrum of the system, which becomes 
ladder-like with the level spacing This structure is 
revealed in photoluminescenc o 59 i 60 and photoconductiv- 
itjiSS spectra as a series of equally spaced peaks. Disorder 
broadens the peaks and makes them unresolved. Below 
we derive a relationship between magnitudes of disorder 
a and bias U, which governs the occurrence of the WSL 
in disordered systems. 

Briefly, our reasoning is as follows (a detailed study 
will be published elsewhere^) . According to the ex- 
change narrowing co ncep t (see the discussion in Sec. IIIII 
and Refs. Ir39ll7lll72ll73ll74li . a quasiparticle confined within 
a chain segment of size Lb sees a disorder of a reduced 
magnitude a j 'y 'Lb, with a the strength of the bare disor- 
der. Therefore, the resulting inhomogeneous broadening 
of the WSL levels can be estimated, similarly to Eq. I|5b|l . 
as 

2(7 /— 

<tb = -7= = <tVTJ . (6) 
V Lb 

In order to resolve the ladder structure, the inhomoge- 
neous width <jb must be smaller than the level spacing 
U. This brings us to the condition, 



which governs the occurrence of the WSL in ID disor- 
dered systems. We will refer to the relationship JJJ) as to 
the limit of strong bias. 

Now, we turn to discussing the opposite sign of in- 
equality (jTJ), which we will name the low or/and moder- 
ate bias limit. At non zero bias, each site gets an addi- 
tional energy —U(n— N/2). Thus, the energy difference 
between the edge sites (magnitude of the total bias) is 
about UN. It is to be compared first of all to the ab- 
sorption line width a* in the absence of bias. Apparently, 
at UN <C cr* , i.e., when the potential drop across the 
whole system is much smaller than the width, the effect 
of bias on the absorption is negligible. On the contrary, 
the bias is expected to broaden the absorption spectrum 
when UN ^S> cr* ■ Indeed, optically dominant (bell- like) 
states in this case will be distributed from [—UN/2 to 
UN/2], giving rise to almost constant absorption within 
this energy range. 

Such a scenario, however, holds as long as UN* < cr*, 
i.e., provided the typical potential drop across the poten- 
tial wells supporting bell-like states is smaller than the 
typical well depth a* in bias-free systems. At UN* > cr* , 
the picture of localization in the potential wells does not 
work anymore. The localization of states now is governed 
by a complicated interplay of disorder and bias, which is 
hard to handle qualitatively. Nevertheless, we can still 
claim that the FWHM of the absorption spectrum will 
be on the order of magnitude of the total bias, UN. On 
further increasing the bias, we fall in the regime of the 
WSL, y/U > a. Our numerical simulation confirm this 
qualitative picture, as shown below. 

Note that the above reasonings can be applied without 
any remarks to both uncorrelated and weakly correlated 
disorder (a < 1). In the limit of strong correlations (a ^S> 
1) the picture requires corrections which we discuss in 
Sec. ESI 




FIG. 2: Absorption spectra of biased chains (U = 0.01) with 
N — 100 sites calculated for two values of the correlation 
exponent a, shown in the plot. The magnitude of disorder 
is cr = 1. Each curve were obtained after averaging over 10 6 
realizations of disorder. 



a < VU , 



(7) 
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B. Low and moderate bias 

In Fig. |21 we plotted the absorption spectra calculated 
for disordered biased chains of N = 100 sites, choos- 
ing the disorder strength a = 1 and the bias magnitude 
U = 0.01. Averaging over 10 6 realizations of disorder 
were performed in Eq. Q). Two values of the correla- 
tion exponents were considered, a = 1 and a = 4. As 
u VU, we are not in the WSL regime. Furthermore, 
at U = 0.01 the overall potential drop UN — 1 is smaller 
than the FWHM in the absence of bias, a* as 2.2 (see 
Fig. As a consequence, the effect of bias is weak, 
leading only to a unnoticeable broadening of the spec- 
tra and smoothing the shape of the doublet (at a = 4) as 
compared to the bias-free conditions (compare to Fig.Q. 
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FIG. 4: Same as in Fig. El but for a bias U = 0.5. 
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FIG. 3: Absorption spectra of biased chains (U = 0.1) with 
N = 100 sites calculated for two values of the correlation 
exponent a, shown in the plot. The magnitude of disorder 
is <r = 1. Each curve were obtained after averaging over 10 6 
realizations of disorder. Insets show enlarged views of the 
spectra within the grey boxes. 

Figure [21 presents the results of the simulations for a 
larger magnitude of the bias, U — 0.1, keeping all other 
parameters unchanged. Still, a > y/U that is not in favor 
of the WSL. However, the total bias UN — 10 is now 
larger than the bias-free FWHM w 2.2. Therefore, for 
both values of the correlation exponent a = 1 and a = 4 
the absorption spectrum shows a plateau-like shape and 
a large broadening, with the FWHM w 10 equal to the 
total bias. These trends are in full agreement with our 
qualitative reasoning presented in the preceding section. 



C. Strong bias 

Aiming to hnd fingerprints of the WSL in the optical 
absorption spectra, we further increased the magnitude 
of the bias. Figure 0] shows the spectra calculated for the 
bias magnitude U — 0.5. As before, chains of N = 100 



FIG. 5: The oscillator strength distribution for the set of 
parameters as in lower panel of Fig. 



sites were used in the simulations and two values of the 
correlation exponent a were considered, a — 1 and a = 4. 
The disorder strength was set to a = 1, and averaging 
over 10 6 realization of disorder was performed in Eq. J3J. 

We observe (the upper panel) that at a = 1 < a c = 2 
the spectrum remains structureless that is consistent 
with the estimate J7J: still a > \/U, i.e., the WSL can 
not be resolved. However, for strong correlations in dis- 
order, when a = 4 > a c = 2, the spectrum presents a 
periodic pattern which is not masked by the stochastic 
disorder fluctuations (see the inset in the lower panel). 
Most important, the period of the modulation is exactly 
equal to U = 0.5. To illustrate more our statement, we 
depicted in Fig.[5]the oscillator strength F v as a function 
of the eigenenergy E v for a = 4. It is clearly seen the pe- 
riodicity in energy with a period equal to U — 0.5 as well 
as the similarity of the corresponding oscillator strengths. 
We checked out that the pattern with the same spacing 
also holds for larger systems, i.e, its period is not a con- 
sequence of finite size effects. Therefore, we claim that 
the periodic pattern found in the simulations results from 
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FIG. 6: Same as in Fig. but for a disorder magnitude a — 
0.2. The inset shows an enlarged view of the center of the 
band. 

the occurrence of the WSL in the energy spectrum of the 
system. 

An explanation of the occurrence of the WSL in the 
strong correlation regime (a = 4) is based on the already 
mentioned in Sec. IIVI fact: the site potential 10) for a 
given realization of disorder is cosine-like (harmonic with 
k = 1), perturbed by a colored noise (harmonics with 
k > 2). This implies that the magnitude of the actual 
disorder is at least by a factor of 2~ Q / 2 smaller than the 
bare value er; the inequality Q turns out to be fulfilled for 
the reduced disorder. In other words, strong correlations 
in disorder facilitate the occurrence of the WSL. 

Finally, we calculated the absorption spectrum profile 
for magnitudes of disorder a — 0.2 and bias U = 0.1, 
when the inequality J7J holds even in the absence of cor- 
relations in disorder (a = 0). The results, obtained for 
a = and a = 5 are depicted in Fig. (upper and lower 
panels, respectively). One observes that the absorption 
spectrum shows a resolved structure for both values of 
a, which, however, is heterogeneous. The central part of 
the spectrum is of our primary interest because it reveals 
a periodic pattern. 

At a = (uncorrelated disorder), the pattern consists 
of equally spaced single peaks with the spacing exactly 
equal to U = 0.1 (see the inset). This again allows us 
to associate the peaks with the occurrence of WSL in 
the energy spectrum of the system. The FWHM of a 
single peak is about 0.06 that is in good agreement with 
Eq. 10. We checked also out that the pattern did not ap- 
pear if a = 0.4 w \/U = 0.45, in full accordance with our 
reasonings presented in Sec. IV Al These results corrob- 
orate those found in Ref. |6^ for random Kronig-Penney 
models. 

In the case of strongly correlated disorder (a = 5), each 
single peak of the central pattern splits into a doublet, as 
is seen from the inset in the lower panel. The origin of the 



doublets can be traced back to the behavior of the un- 
biased system, where the absorption line shape exhibits 
a doublet structure provided a > a c = 2 (see Sec. IIVI 
and Ref. I66T) . Figure HO points out that the splitting also 
occurs in the presence of bias. 

To conclude we comment on the peculiarities of the red 
and blue sides of the absorption spectrum. We associate 
them with finite size efects. Indeed, these parts of the 
energy spectrum are formed by the states localized close 
to the system ends. Because of that, the corresponding 
eigenfunctions differs from those at the band center. 



VI. SUMMARY AND CONCLUDING 
REMARKS 

We studied numerically the linear optical absorption 
of a quasiparticle moving on a ID disordered lattice sub- 
jected to a linear bias of magnitude U. The random site 
potential was set to have a power-law spectral density 
S(k) ~ l/k a , which gives rise to long-range correlations 
in site energies. 

The absorption spectrum of the unbiased lattice (U = 
0) was found to present a single peak in the weakly- 
correlated limit, when the correlation exponent a < a c = 
2, with a c — 2 the critical value for the occurrence of 
a phase of extended states in the center of the band. 
For strongly-correlated disorder, at a > a c — 2, the ab- 
sorption lineshape turned out to exhibit a double-peaked 
structure, characteristic for this type of long-range cor- 
relations*^ 

Switching the bias on does not change much the out- 
lined behavior provided the magnitude of overall bias U N 
does not exceed the absorption bandwidth of the unbi- 
ased system. At higher magnitudes of bias, the absorp- 
tion spectrum starts to broaden, independently of the 
magnitude of the correlation exponent a. Its profile gets 
a top hat shape, being almost flat at the center of the 
absorption band and having the FWHM on the order of 
UN. Such a scenario holds as long as the disorder mag- 
nitude a < VU. 

On further increasing the magnitude of the bias, a pe- 
riodic pattern is found to build up at the center of the 
(already wide) absorption band. Its period is equal to U, 
as for the Wanier-Stark ladder in an ideal lattice, and in- 
dependent of the system size N . Therefore, we attribute 
the pattern found to the Wannier-Stark quantization of 
the energy spectrum in the disordered lattice. 

The occurrence of the Wannier-Stark pattern is facili- 
tated by the presence of correlations in disorder. In the 
limit of strong correlations (a > a c = 2), each Wannier- 
Stark level represents a doublet, reflecting the doublet 
structure of the absorption spectrum of the unbiased sys- 
tem. 

To conclude we note that in our study we did not 
explicitly relate the bias to the presence of an external 
uniform electric field. Therefore, our conclusions can be 
equally applied both to ID disordered electrons, moving 
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in a uniform electric field, and to ID disordered Frenkel 
excitons, where energetic bias can be an intrinsic prop- 
erty of the system. Dendrimers represent one of the ex- 
amples^ 
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